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PACS numbers: 13.60Le, 13.85Qk 



Since it was first observed, charmonium (cc) produc- 
tion in hadronic collisions has been a subject of consid- 
erable theoretical interest. Recent approaches to under- 
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standing charmonium production make use of nonrela- 
tivistic QCD [J, Q to calculate hadro-production rates 
at the Tevatron and elsewhere. While most experimental 
observations of charmonium production consist of J/tp 
measurements, a significant contribution of J/tp produc- 
tion is indirect, resulting from the decay of higher mass 
states 3[. In particular, the radiative decay of the Xcj 
states [4J accounts for a significant fraction of the J/tp 
production seen in hadronic collisions, and any calcula- 
tion of J/tp production must include Xcj production as 
well. 

Measurements of hadronic \ c j production have been 
made in a variety of beam types and energies @ by ob- 
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serving the decay process XcJ ~ * J/ip 7- Experimental 
results available until now have suffered from large sta- 
tistical uncertainties, and no measurement has had the 
precision to test the consistency of the cross-section ra- 
tio cr Xc2 /<7y cl with the simple spin-state counting expec- 
tation of | for XcJ mesons that are directly produced 
in the interaction 0. Knowledge of this ratio is needed 
in calculations of J/ip production through radiative XcJ 
decay, and can be an important standard for comparing 
production models. 

In this letter, we report a measurement of the rela- 
tive cross section times branching fractions of the Xd 
and Xc2 mesons produced in pp collisions at a center of 
mass energy of 1.96 TeV using the CDF II detector at 
the Fermilab Tevatron. We study the inclusive process 
PP -> XcjX, where XcJ -> J/ipl, and J/ip -> fi + fj,~, in 
a data sample with a time-integrated luminosity of 1.1 
fb _1 . The final state photon is reconstructed through 
its conversion into e + e~, which provides the mass res- 
olution needed to distinguish the Xd and Xci states. 
The spatial resolution of the [i + (J.~ vertex allows sepa- 
ration of prompt Xd production from events where the 
XcJ meson is a i?-hadron decay product. We measure 
the ratio of the cross section times branching fraction 

R P = c r Xc2 K(xc2 -> J/^l)/v Xcl B(xci -> J/ipy) for 
promptly produced XcJ mesons. In addition, we obtain 
a result for the analogous quantity in B decay events, 
R B = o B B(B -> X c2X)B(xc2 -> JM>i)I<t b B(B 
Xci X)B(xci - > J/ipj)i which provides a measurement 
of B(B -> X c2 X)/B(B -> Xd X) for the B hadrons pro- 
duced in the Tevatron environment. 

This analysis makes use of the tracking, muon identi- 
fication, and trigger systems. The CDF II detector has 
been described in detail elsewhere [1, [§]. The tracking 
system consists of a seven-layer silicon microstrip detec- 
tor and an open-cell drift chamber (COT) that operate 
inside a solenoid with a 1.4 T magnetic field. Muon 
candidates from the decay J/ip — > are identi- 

fied by two sets of drift chambers located outside the 
electromagnetic and hadronic calorimeters. The central 
muon chambers cover the pseudorapidity region 1 77 1 < 0.6, 
and are sensitive to muons having transverse momentum 
p T > 1.4 GeV/c [l(|. A second muon system covers the 
region 0.6 < |7y| < 1.0 and is sensitive to muons having 
Pt > 2.0 GeV/c. Muon triggering and identification are 
based on matching tracks measured in the muon system 
to COT tracks. 

The analysis of the data begins with a selection of well 
measured J/ip — > candidates. These are selected 

by requiring events that contain two oppositely charged 
muon candidates, each with a match between the COT 
and muon chamber tracks. We also require that both 
muon tracks have measurements in at least three layers 
of the silicon detector and a two-track invariant mass 
within ±80 MeV/c 2 of the world-average J/ip mass Q. 
The J/ip candidates are required to fall within kinematic 
bounds of p T (J/ip) > 4.0 GeV/c and \r)(J/ip)\ < 1.0, 
which correspond to the approximate limits of our ac- 



ceptance. A simultaneous mass and vertex constrained 
fit is performed on the muon tracks, where the dimuon 
mass is constrained to the world-average J/ip mass. 

The search for photon conversion candidates begins 
with a scan of all additional tracks with pt > 400 MeV/c 
found in each J/ip event. Two oppositely charged tracks 
are each assigned the electron mass, and have their track 
parameters recalculated by subjecting them and their un- 
certainties to a fit that has constraints consistent with 
the photon conversion hypothesis. Specifically, the two 
tracks are constrained to be parallel at their point of in- 
tersection, and the momentum vector of the pair is con- 
strained to originate from the dimuon vertex. A dis- 
placement of 12.0 cm or more from the beam line in the 
direction of the track pair's transverse momentum is re- 
quired to omit conversions whose momentum is poorly 
measured due to bremsstrahlung in the inner detector 
material. We also require Pt(j) > 1-0 GeV/c. Finally, 
a constrained fit is performed on the four tracks that 
combines the J/ip mass constraint with the photon con- 
version hypothesis. The invariant mass distribution of 
all J/ijj"f combinations is shown in Fig. [TJ which clearly 
demonstrates that the J/ip 7 mass resolution acheived by 
this technique is sufficient to resolve the XcJ states. 

The lifetime of B hadrons allows the transverse dis- 
placement of the dimuon vertex from the beamline to be 
used as a tool for their identification. Since any J/ip"f 
combination that originates from B decay represents only 
a partial reconstruction of the B hadron, the proper life- 
time is not directly measurable. We therefore use the 
quantity ct = Lx M(J/ip) x F(pr(J/ip))/pr(J/ip) where 
M(J/ip) and pr{J/ip) are the mass and transverse mo- 
mentum, respectively, of the J/ip candidate, and L is the 
measured displacement of the dimuon vertex in the di- 
rection of pt{J /ip)- The correction factor, F(pr(J/ip)), 
was obtained by a Monte Carlo simulation of B hadron 
decay [3| , and provides an average correction between the 
measured displacement and the lifetime of the decaying 
B hadron. 

Mass resolution, acceptance, and reconstruction effi- 
ciencies for J/ip 7 final states of different invariant mass 
have been studied with a Monte Carlo simulation that 
generates events uniformly in rapidity and with a trans- 
verse momentum distribution that matches the measured 
distribution for J/ip events The simulated events 
were processed through our reconstruction and analysis 
algorithms, and provided templates for the expected sig- 
nal shape of the final J/ip 7 invariant mass distribution 
as a function of pr( J/ ip)- In particular, the simulated 
events enabled an estimate of the mass resolution and 
radiative tail due to scattering and radiation by the con- 
version electrons as they passed through the material in 
the detector. 

We used an unbinncd likelihood fit to calculate the 
yield of XcJ events for both prompt and -B-decay pro- 
duction processes. The probability density function used 
for the fit is a function of both J/ip 7 invariant mass 
and ct. Independent signal and background distribu- 
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tions are calculated for both processes. The mass dis- 
tributions of the signals are constrained to the templates 
obtained through simulation. The mass distributions of 
the backgrounds are modeled by polynomials, and the 
probability density function for each event uses the cal- 
culated uncertainty on the invariant mass and ct. The 
ct distribution is used to separate the production pro- 
cesses, and is modeled as a sum of prompt (Gaussian 
resolution) and S-decay (exponential, convoluted with 
resolution) contributions. Our fit to the data gives an 
event yield (N XcJ ) of iV Xc0 = 41 ± 20, N Xel = 2143 ± 60, 
and N Xa2 = 1035 ± 40 for promptly produced events. 
For B decay events, the yields are N Xc0 — 29 ± 16, 
N Xcl = 384 ± 35, and N Xa2 = 66 ± 16. Projections of 
the unbinned likelihood fit are overlaid onto the mass 
and ct distributions shown in Figs. [1] and [51 The rela- 
tively small yield of \c0 candidates is due to the small 
branching fraction into the J/ip 7 final state Q , and is 
the reason they were not used in the subsequent analysis. 

For our acceptance calculation, we have analyzed our 
simulated events assuming every photon converted, and 
created electron-positron pairs according to the Bethe- 
Heitler distribution [TJ. We then counted the number 
of events that would have been accepted if all final state 
products were to pass our kinematic requirements. Our 
low transverse momentum limit of 400 MeV/c for the 
electrons results in a dependence of our acceptance on 
the invariant mass of the parent. The overall ratio of the 
Xci and Xc2 meson acceptances (e Xc2 / e Xcl ) is listed in Ta- 
ble U for several ranges of Pt(J/iP)- The acceptance ratio 
is then combined with our yield ratios to provide mea- 
surements of R p and Rb for several ranges of pt{J/^)- 

TABLE I: The acceptance ratio and ratios of cross section 
times branching fractions of the Xcj states for the prompt 
events and B decay events. Uncertainties listed are statistical 
only. 



p T {JM) (GeV/c) 


e xc2/ e xd Rp R B 


4-6 


1.27 ± 0.01 0.457±0.039 0.150±0.087 


6-8 


1.17 ±0.01 0.384±0.034 0.080±0.094 


8-10 


1.14 ± 0.01 0.455±0.053 0.116±0.070 


> 10 


1.10 ±0.01 0.309±0.045 0.197±0.082 



> 4 1.23 ± 0.01 0.395±0.016 0.143±0.042 



tribution given by /(6> w ) oc 1 + a cos 2 where 9^ 
is the polar angle of the /i + in the J/tp rest frame, and 
we take a = 0.13 ± 0.15 as was done for a separate J/tp 
cross-section measurement Q. A variation of a by this 
uncertainty for one of the two XcJ states implies a relative 
shift of ±0.030 on our reconstruction efficiency ratio. 

We have also considered the sources of systematic un- 
certainty in the yield ratio calculation. We varied the 
invariant mass signal shape obtained from the simulation 
within the uncertainty of its parametrization, and found 
that the relative shift in the yield ratio is ±0.005. The 
uncertainty on parameters used in the ct definition corre- 
sponds to a variation of the B fraction of the J/ipj events 
of ±0.007 of its value, giving a systematic uncertainty of 
±0.002(±0.010) in the yield ratio for the prompt(B de- 
cay) sample. Finally, we explored the possibility that 
our data contain partially reconstructed h c — > J/t[)ir°, 
7T° — ► 77 events. This has been studied by simulating 
this process, parametrizing the resulting J/tjj 7 invariant 
mass distribution, and including this possible background 
in our likelihood fits for the signal yield. The possible h c 
background contribution was found to be negligible in 
this data sample, so no systematic uncertainty was as- 
signed for this process. 

Differences in the two states due to production angu- 
lar or px distributions would require different production 
mechanisms for the Xd an d Xc2 mesons, and are, there- 
fore, considered to be unlikely. Consequently, we did not 
assign a systematic uncertainty on the XcJ acceptance 
due to production dynamics. A summary of the system- 
atic uncertainties on the cross section times branching 
fraction ratio is listed in Table [TTJ The individual un- 
certainties are combined in quadrature to give the total 
systematic uncertainty. 

TABLE II: Relative systematic uncertainties on R p and Rb- 



Effect Uncertainty 

Simulation Sample Size ±0.005 

Photon Conversion Simulation ±0.020 

Polarization Effects ±0.030 

Invariant Mass Resolution ±0.005 

Prompt/B Separation ±0.002(±0.010 for B) 

Total ±0.037 (±0.038 for~B) 



Several systematic effects that might change the recon- 
struction efficiency ratio e X c2/ e xci were studied. First, 
the simulated event sample size used for the acceptance 
calculation provides an overall relative uncertainty of 
±0.005 on the ratio. A comparison between the full 
event simulation/reconstruction and the simpler simu- 
lation based on the electron energy distribution yields a 
relative systematic uncertainty of ±0.020. Another ef- 
fect considered is that polarization of one of the two XcJ 
states would also introduce a systematic shift. We have 
evaluated the effect of having one state decay with a dis- 



Our final result on the relative rate of production 
for promptly produced XcJ states is R p = 0.395 ± 
0.016(atat.) ± 0M5(sys.) for X c.j with pr(J/ip) > 4 
GeV/c and pt{i) > 1 GeV/c. For XcJ resulting from 
B decay over the same kinematic range we find Rb — 
0.143 ± 0M2(stat.) ± 0.005(sys.). These results provide 
the most precise measurement of the XcJ production ratio 
obtained in any hadronic interactions. Conversion of this 
measurement into the direct cross-section ratio a Xc2 /<J Xc i 
requires a knowledge of the branching fractions, which 
are not measured in this experiment, and a small correc- 
tion due to ip{2S) — > Xcj 7 decays. Based on the existing 
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measurement of the prompt ip(2S) cross section 12j, and 
the Xcj contribution to the prompt J/ip production cross 
section Q, we estimate that 4.0 ± 1.0(5.0 ± 1.0)% of our 
prompt Xci{Xc2) sample is due to decay of promptly pro- 
duced ip(2S) mesons. 

Prior measurements of the prompt cross-section ratio 
have been severely limited in their precision due to the 
statistical uncertainties inherent to small data samples 
The relative precision of previous measurements has 
typically been approximately 30% on the cross-section 
ratio, and provides weak guidance for production models. 
This work, combined with the best branching-fraction 



ratio measurement R 



B( Xc i ^ J/V7)/B(Xc2 



J/ipi) = 1.91 ±0.10 available [13|, gives Rp Rj/ipj — 
0.75±0M(stat.)±0M(sys.)±0.04:(B.F.), where the last 
term in the uncertainty is due to the branching-fraction 
ratio uncertainty. This level of precision should serve 
to inform any future developments in the calculation of 
hadronic charmonium production. 
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FIG. 1: The J/ip"f mass distribution (points) with the pro- 
jection of the likelihood fit overlaid on the data. The masses 
of the Xcj mesons and the contributions of the signal and 
background components are indicated. 
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FIG. 2: The ct distribution (points) for events in the \ a \ 
(a) and Xc2 (b) mass ranges. The projection of the fit is 
overlaid on the data, with the contribution of each signal and 
background component indicated. 



